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Abstract 

For a 4-dimensional 3-parametric toy Hamiltonian H(a, b, c) we construct the domain 
T> of couplings in which the eigenvalues E n remain real (i.e., in principle, observable). 
A relationship is found between the reflection symmetry of the spectrum (i.e., its 
Dunne's and Shifman's self-duality Ej = —E^ + i_j at iV = 4) and a geometric 
symmetry of the physical domain T>. Simultaneously, both remain unbroken at a = c. 
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1 Introduction 



Reflection symmetry 

E j = -E N + 1 .. i , J = 1,2,... (1) 

between the low- and high-excitation parts of the energy spectra is one of features 
of all the many-body Hamiltonians which are quadratic in the creation and anni- 
hilation operators pQ. In a different physical context of quasi-exactly solvable local 
potentials, Dunne and Shifman [2] found the same form of symmetry and called it 
self-duality of the spectrum. Recently, it was rather surprising for us to detect the 
emergence of the same spectral symmetry also during the study of the strong- 
coupling version of a family of finite-dimensional phenomenological non-Hermitian 
(often called VT— symmetric) chain- model Hamiltonians [3]. 

The latter study revealed that the additional symmetry (pQ) appeared in connection 
with a thorough simplification of the spectra and also with a facilitated feasibility 
of an algebraic determination of the positions of certain maximal physical couplings. 
Formally, the latter maximal couplings were defined as vertices of the boundary 
dT> of the domain D where the free parameters of our models remained physically 
acceptable, i.e., compatible with the reality and observability of the energies. As 
long as the domain D proved bounded, the observability of the related energies can 
be called "conditional" [4J. In our subsequent deeper study [5] of the same class 
of models with symmetry (pQ) we succeeded in a further constructive extension of 
the previous result and in a clarification of the geometric structure of the horizons 
dT> in non-empty vicinities of all the above-mentioned vertices. In our present brief 
communication we feel inspired by the latter experience. 

Our purpose is to throw new light on the peculiar apparent correspondence be- 
tween the spectral symmetry flTJ and the geometric shape of the boundary &D. An 
explicit description of the latter manifold will be performed here for a specific, modi- 
fied model where the up-down symmetry ([1]) will be manifestly broken. We shall see, 
in essence, that after the violation of symmetry ([1]) also the mysterious simplicity of 
the shape of the surface dT> gets lost. For the sake of clarity of our argument, we 
shall only consider the construction of &D for the exactly solvable three-parametric 
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The problem differs from its two-parametric b = c predecessor of ref. [3] (cf. also a 
few related marginal remarks in [6]) by a non- vanishing measure / > of the extent 
of the asymmetry which enters the upper coupling defined as b = b(c, f) = y 'c 1 + p. 

We shall not analyze the possibilities of a transition to the higher- dimensional 
models here. A thorough discussion of such a perspective has already been presented 
in both our preceding papers on self-dual cases [31 [5] . Basically, we saw there that 
the key information about the structure of the spectra is already provided by the 
models of the lowest dimensions. The growth of the matrix dimension appeared to 
provide merely technical alterations of the overall picture. 

On this background, we now intend to pay more attention to the most elementary 
up-down-asymmetric model ([2]). We shall apply computer-assisted linear algebra 
and formulate some of its consequences in sections [2] and [3J with a brief summary in 
section HI 



2 Physical domain T> (graphical approach) 

Let us start by a terminological remark concerning eq. ([2]) and all its generalizations 
with parameters inside D. Its reason is that the accepted denotation of similar 
models characterized by a non-standard transposition asymmetry differs in different 
applications. We will prefer calling our H quasi-Hermitian (a term typical for nuclear 
physics (Tj). Still, several authors also use the equivalent names of pseudo-Hermitian 
H (predominantly in the context of linear algebra and mathematics [8J) or crypto- 
Hermitian H (this nice and most self-explanatory concept appeared recently in the 
context of gauge models [9]). Last but not least, another, less strict nickname of 
VT— symmetric H as coined by Carl Bender [10] became most popular as particularly 
appealing in field theory. 
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2.1 The allowed range of A(a,c,f). 

Returning now to our particular model let us recollect that its energies coincide with 
the roots of its secular determinant, 

det(H-EI)=Q, (3) 

where we decided to set b 2 = c 2 + f 2 , with a new measure of the up-down asymmetry 
/ > replacing the original coupling parameter b. We may express eq. ([3]) in its 
explicit polynomial form 

^2,3,4 - A{a, c, /) E 2 m + C(a, c, f) = 4 f 2 E w (4) 

where we have A(a, c, f) = 10 — a 2 — 2 c 2 — f 2 while 

C(a, cJ) = 9 + Qc 2 -9a 2 + 3f 2 + c 4 + fc 2 

or, in compact form, 

C(a, c, /) = (3 + c 2 ) 2 + f (3 + c 2 ) - 9 a 2 = (3 + c 2 + / 2 /2) 2 - (9 a 2 + f 4 /A) = 

= (8 - A/2 - a 2 /2) 2 - (9 a 2 + f/A) . 
We observe that at any positive / > and negative (or vanishing) A < 0, the left- 
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Figure 1: Graphical solution of our secular eq. (jUat A > 0). 

hand-side curve of eq. (jlj) can only have at most two real intersections = Ei^(f) 
with the right-hand-side straight line. This means that the spectrum always contains 
at least two complex conjugate energy roots and we stay off V. 
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If we allow A < and vanishing / = 0, the two complex roots become real only 
in the limit A — > 0. Thus, the point of the boundary &D of the physical domain of 
parameters would be achieved, in full agreement with our older self-dual results [3]. 
Now, we decided to assume that the self-duality is manifestly broken, / > 0, so that 
we may conclude that we have always to consider just the strictly positive values of 
the polynomial 

A(a,c,f)>0, {a,c,f)eV. (5) 

Under this assumption, the graphical solution of eq. (TjJ is sampled, in arbitrary units, 
in Figure [H At a variable / > we shall further require that inside £>, the left-hand 
side superposition of the two (viz., quadratic and quartic) parabolic curves of eq. (jl]) 
has strictly four real intersections with the right-hand side straight line. From this 
assumption one immediately deduces that we must have a negative E±(f) > Ei(0) 
and E 2 {f) < E 2 (0) and a positive E 3 {f) < E 3 {0) and E 4 {f) > E 4 {0). Hence, the 
two rightmost roots E 3 and E 4 can never merge and, subsequently, they cannot form 
a complex conjugate doublet unless / = A — 0. In terms of the geometry of T>, the 
transition to / 7^ immediately destroys the multidimensional hedge-hog-shaped 
form of the surface dT> as observed at / = 

2.2 The allowed range of C(a,c,f). 

An inspection of Figure [T] reveals that at E = 0, one of the curves acquires the value 
of C(a,c,f) which may be, in principle, positive, vanishing or negative. Once we 
keep the quantity / = /o > fixed, the variation of a and/or c just changes the val- 
ues of the functions C(a, c, / ) and A(a, c, / ). In the two-parametric space of a and 
c we may move along a line of a constant A(a, c, /o) and observe that any decrease 
of C(a,c,fa) leads to the decrease of the distance between the intersections E 2 (f ) 
and E 3 (f Q ). In principle, these two values will merge and subsequently complexify 
below certain f — dependent value C( miTOts )(a, c, / ) < 0. Similarly, the growth of 
the value of the function C(a, c, / ) will imply the decrease of the distance between 
the two leftmost energies Ei(f ) and E 2 (fo)- The latter pair will merge (and subse- 
quently complexify) at an upper bound C(p/„ s )(o, c, fo) which, by the way, need not 
be necessarily positive. We may conclude that all the four energies remain real iff 

C (minus) (a, c, /) < C(a, c, /) < C(pi us ) (a, c, /) , (a, c, /) G V . (6) 
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Both these bounds follow also algebraically from our secular equation (j3J) of course. 
Their occurrence corresponds to the confluence of the above-mentioned pairs of en- 
ergies at a point of the boundary &D of T>. 

In the graphical language this means that inside all the interior of T>, the graph of 
the secular polynomial Y(E) = det(H - E) = E 4 - A E 2 - 4/ 2 E + C will intersect 
the real line four times, achieving its leftmost, negative minimum at a negative real 
variable Z( min ) G (Ei,E 2 ), 

Y{z(min)) = z\min) ~ A Z( m in) ~ Z( m in) + C( p i us ) < . (7) 

The subsequent positive maximum of Y(z) will occur at a larger Z( max ) G (E%, E 3 ), 

Y(z( max )) = Z( max j — A Z 2 ^^ — Af 2 Zf^nax) + C( m i nus ) > . (8) 

As long as we showed above that the two rightmost energy roots E3 and E4 can 
never merge, the remaining, third extreme of the function Y(z) at some z > z max 
may and will be ignored here as redundant. 
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Figure 2: Graphical solution of eq. (jHJ) at A > 0. 

With the variation of a, c and /, we shall always have Y(0) = C while the quan- 
tities Z( m i n ) < z^ax) correspond to the respective leftmost minimum and subsequent 
maximum of the difference curve Y(z) (in Figured! you may imagine that / varies). 
Algebraically, the two quantities Z( m j n ) and Z( max ) coincide with the two real roots of 
the third-order polynomial dY/d z , i.e., of the cubic equation 

\min J \ ™ J 
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(cf. Figure [2]). Once we know these auxiliary values, we can define the required 
bounds from their definition 



C f . ,=Afz f \+Az 2 f \ (10) 

/ minus \ J l max \ I m.n.rr. \ I m.ncr. \ > f 



if. \ = <i J Zf \ -f AZ f \ 

I minus \ <> I max \ I max \ 

\plus J \min J I min J 



which, after an algebraic simplification using eq. (jH]), degenerates to the simpler 
formula 

= ,+3f Z f \ . (II 



'/ ■ \ = -2/ \ + of 2/ \ 

/ minus \ O mm J I max\ 

\plus J ^ [mini [mini 



We may conclude that since the only negative component of the polynomial C(a, c, /) 
is —9 a 2 , the latter bounds in effect fix the allowed range of the parameter a for any 
given pair of c and /. In the other words, the middle energy levels E 2 and E 3 
can merge iff the value of a proves sufficiently large. Moreover, since A > 0, the 
quadruple merger of E\ and E 2 and E$ with E± would require a return to / = and 
proves entirely excluded at / ^ 0. 



3 Physical domain T> (analytic approach) 

In our previous studies [31 [5] of iV— dimensional chain models it has been conjec- 
tured and verified that at all N, the postulate of a self-duality of the spectrum 
Ej = —En + \^j could perceivably simplify the construction as well as geometric 
characterization of the shape of the manifold T>. In the preceding section [2] an 
independent constructive support of the plausibility of such a relationship has been 
given via our three-parameter non-self dual example ([2]) . In what follows we intend to 
complement these results by paying attention to a simplification and non-numerical 
description of the boundary dT>. This might prove relevant, e.g., in the quantita- 
tive context of the strong-coupling and/or singular forms of perturbation theory [5] 
as well as in a more qualitative setting of a possible quantum version of theory of 
catastrophes [1]. 

3.1 The allowed range of the asymmetry /. 

We assumed the knowledge of the input energy-degeneracy-determining auxiliary 
quantities Z( m j n ) and Z( max ) which can both be defined by the Cardano formulae in 
principle. Fortunately, a perceivably less painful recipe can be recommended. In 
a preparatory step, an inspection of eq. (cf. also Figure [2]) reveals that the 
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pair of our auxiliary quantities Z( m i n ) < Z( max ) remains real iff / € (0, f (upper))- 
Here, the numerical value of the upper bound f( U pper) annihilates the polynomial 
Z(z) = 4 z 3 — 2 A z — 4 f 2 at its left maximum where dZj d z = at z = Z( upper y This 
means that 

2 [A [A 3 \ 1/4 

12 Z(u VV er) = 2 v4 , Z(u pper ) = — W — , f (upper) = I I 

Thus, we are allowed to re-parametrize / — > f in such a way that 

/ 2 = /(Uer) COS V > V = <f( A , f) e (°> */ 2 ) . 

Moreover, using the direct insertion in eq. (Q we easily verify that its closed solutions 
simply read 

2A /7T - 



Z(min) V Q 




3 V 3 
2^4 /^ + tt 




*^(ma:r) ^'3 i g 

Both of them are negative and they only have to satisfy the pair of constraints 

Z(min) = Z( min) (A,f) G (-y/A/2, — \J A/6) , (12) 
Z(max) = Z(rnax)(A, if) G (-\J A/6, 0) . (13) 

This enables us to conclude that all the above-mentioned open interval of f may 
be treated as lying inside T>. The angle if is, therefore, a better measure of the 
breakdown of the self-duality of our toy Hamiltonian H . Next, let us show that and 
why also its further coupling constants should be reparametrized. 

3.2 A reparametrization (a, c, /) — ► (a, 5, if) of H and V 

Let us replace A — ► a in such a way that A — 10 sin 2 a with a G (0, vr/2). Moreover, 
let us also reparametrize a — ► 5 with a 2 = 10 cos 2 a sin 2 5 and 5 G (0, 7r/2). This 
implies that the resulting formula 2 c 2 + / 2 = 10 cos 2 a cos 2 5 may be inserted in the 
definition of 

C(a, c, f) = (3 + 5 cos 2 a cos 2 5) 2 - 90 cos 2 a sin 2 5 - f A /A = 
= (l2 + 5 cos 2 a cos 2 5) 2 - 90 cos 2 a - 135 - / 4 / 4 • 
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In parallel, we may also reparametrize eq. (TTTj) . 

r /7T + ,0\ 1 

(14) 



7T ± (f' 

cos I — - — I — COS (f 



100 . 4 fn± V 
C ( . \ — — — sin a cos — - — 

I minus | Q \ Q 

[plus I \ <J 

which is independent of 5. Thus, we may introduce another pair of shifted and safely 
positive boundaries 

B ( \=C, . x + 90 cos 2 a + 135 + f 4 /A = B ,. \(a,cp) 

| minus \ I minus \ I rn ' irius \ v * 

I plus J \ plus J \ plus J 

and arrive at the final and amazingly compact form of our key constraint (jH]) , 



B(rninus) (a, <f) < 12 + 5 cos 2 a cos 2 5 < y B {plus) (a, ip) , 8eV. (15) 

Our construction of the quasi-Hermiticity [7] domain V is completed since eq. ( !T5|) 
is its definition. In the limit / — > 0, this equation also has been checked to specify 
the simpler domain T> described in paper [3]. 

At / > our present generalized definition (TIB]) of T> is algebraic and rigorous 
but its geometric interpretation is not too transparent. Indeed, although our two 
inequalities specify the allowed, physical range of S for any given pair of parameters 
(a, if), the interval can prove empty. This would mean that the domain T> does 
not contain any point with pre-selected (a, ip). Unfortunately, an account of similar 
subtleties already lies a bit beyond the scope of our present brief communication. 



4 Summary 

Any non-Hermitian VT— symmetric quantum Hamiltonian H remains physical in a 
domain of parameters T> where the spectrum is real and, hence, measurable. In refs. 
[31 EJ we revealed that several properties of a family of tridiagonal matrix Hamil- 
tonian models of this type become exceptionally transparent after an imposition of 
the up-down symmetry or "self-duality" requirement ([1]). In our present study we 
complemented this observation by a test of the consequences of a manifest violation 
of the self-duality. 

We saw that once the real measure f 2 of the violation of self-duality becomes 
different from zero, our first nontrivial, up-down-asymmetric model ([2]) (skipped in 
[6] as complicated) offers a challenging eigenvalue problem. We have shown that its 
discussion and the non-numerical construction of its physical domain of parameters 
T> still remains feasible and compact. 
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We believe that our present construction of a toy model with non-numerically 
tractable physical horizon &D will further improve our understanding of the mathe- 
matics which underlies the so called "conditional-solvability" phenomenon [I] as well 
as the analyses of the so called "quantum catastrophes" [U [11]. At present, the 
•practical role of these fairly fresh mathematical concepts finds new and new physical 
applications, the number of which is rapidly increasing. Besides their above-stressed 
innovative role in quantum theory and particle physics [10], it is worth noting, in the 
conclusion, that their use also currently inspired new progress in the areas as remote 
as relativistic cosmology [8] and phenomenological magnetohydrodynamics [T2] . 
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Figure captions 

Figure 1. Graphical solution of our secular eq. ((4]) at A > 0. 
Figure 2. Graphical solution of eq. ([9]) at A > 0. 
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